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Abstract 

Haematococcus pluvialis cells predominantly remain in the macrozooid stage under favourable environmental condi- 
tions but are rapidly differentiated into haematocysts upon exposure to various environmental stresses. Haematocysts 
are characterized by massive accumulations of astaxanthin sequestered in cytosolic oil globules. Lipidomic analyses 
revealed that synthesis of the storage lipid triacylglycerol (TAG) was substantially stimulated under high irradiance. 
Simultaneously, remodelling of membrane glycerolipids occurred as a result of dramatic reductions in chloroplast 
membrane glycolipids but remained unchanged or declined slightly in extraplastidic membrane glycerolipids. De novo 
assembly of transcriptomes revealed the genomic and metabolic features of this unsequenced microalga. Comparative 
transcriptomic analysis showed that so-called resting cells (haematocysts) may be more active than fast-growing veg- 
etative cells (macrozooids) regarding metabolic pathways and functions. Comparative transcriptomic analyses of asta- 
xanthin biosynthesis suggested that the non-mevalonate pathway mediated the synthesis of isopentenyl diphosphate, 
as the majority of genes involved in subsequent astaxanthin biosynthesis were substantially up-regulated under high 
irradiance, with the genes encoding phytoene synthase, phytoene desaturase, and p-carotene hydroxylase identified 
as the most prominent regulatory components. Accumulation of TAG under high irradiance was attributed to moderate 
up-regulation of de novo fatty acid biosynthesis at the gene level as well as to moderate elevation of the TAG assem- 
bly pathways. Additionally, inferred from transcriptomic differentiation, an increase in reactive oxygen species (ROS) 
scavenging activity, a decrease in ROS production, and the relaxation of over-reduction of the photosynthetic electron 
transport chain will work together to protect against photooxidative stress in H. pluvialis under high irradiance. 

Keywords: Astaxanthin, high irradiance, lipidome, photooxidative stress, transcriptome, triacylglycerol. 



Introduction 

Ever since it was first described by J. von Flotow in 1844, 
Haematococcus pluvialis (Chlorophyceae, order Vovlocales), 
a unicellular freshwater microalga, has drawn significant 



scientific attention due to its potential as a natural resource of 
the high-value carotenoid astaxanthin (Lorenz and Cysewski, 
2000; Cardozo et al. , 2007; Leonardi et al. , 201 1). Astaxanthin 
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is one of the most powerful natural antioxidants and is a 
natural food colourant. Therefore, it is of great commercial 
importance in the nutraceutical and aquaculture industries 
(Lorenz and Cysewski, 2000). In H. pluvialis cells, astaxanthin 
is present in the form of fatty acyl mono- or diesters under 
various environmental stress conditions such as nutrient 
deprivation, high irradiation, and high salt concentrations. 
Upon exposure to these stressors, the green motile flagellated 
and non-motile palmella forms of H. pluvialis cells rapidly 
transform into resting cells called haematocysts or aplano- 
spores. This change is accompanied by massive accumulation 
of astaxanthin in cytosolic lipid bodies, turning the cells a 
distinct red colour (Boussiba, 2000; Grunewald et al, 2001; 
Steinbrenner and Linden, 2003; Damiani et al, 2010). The 
haematocysts can accumulate this versatile secondary carot- 
enoid at up to 4% of cellular dry weight (DW) which is two 
orders of magnitude greater than that found in yeasts and 
higher plants (Boussiba, 2000; Lorenz and Cysewski, 2000; 
Steinbrenner and Linden, 2001; Wang et al, 2009a). 

Astaxanthin accumulation is generally assumed to be 
a survival strategy adopted by some microalgae grown in 
stressful environments, acting as a sunscreen against exces- 
sive irradiance (Hagen et al, 1994), an antioxidant prevent- 
ing production of reactive oxygen species (ROS) (Kobayashi 
et al, 1997; Li et al, 2008), and a carbon or energy reserve 
(Quin et al, 2008) (for a review, see Lemoine and Schoefs, 
2010). As a 40-carbon isoprenoid compound, astaxanthin is 
synthesized from two five-carbon isoprene precursors, iso- 
pentenyl diphosphate (IPP) and dimethylallyl diphosphate 
(DMAPP). Sequential addition of three IPP molecules to 
DMAPP yields 20-carbon geranylgeranyl pyrophosphate 
(GGPP), two molecules of which are condensed head-to-tail 
to tetraterpene phytoene, the precursor for all carotenoids. 
The genes encoding a suite of enzymes that catalyse the for- 
mation of astaxanthin from phytoene have been cloned and 
characterized in H. pluvialis (Steinbrenner and Linden, 2003; 
Steinbrenner and Sandmann, 2006; Vidhyavathi et al, 2008); 
however, the genetic details of the biosynthetic route(s) lead- 
ing to the formation of IPP and DMAPP have yet to be 
elucidated. Two unrelated IPP biosynthesis pathways — the 
mevalonate (MVA) pathway and the 2-C-methyl-D-erythritol 
4-phosphate (MEP) pathway — have been described in vas- 
cular plants. Although it has been suggested that the MEP 
pathway is responsible for carotenoid biosynthesis in vascu- 
lar plants and eukaryotic algae, direct evidence pointing to 
its role in astaxanthin biosynthesis in H. pluvialis is lacking 
at the gene level (Sun et al., 1998; Grunewald et al, 2000; 
Hagen and Grunewald, 2000). A genome-wide survey, along 
with differential gene expression analysis, would be neces- 
sary to identify all genes involved in astaxanthin biosynthesis 
in this organism. A comparison of their expression profiles 
would then be required in order to facilitate identification of 
key regulatory genes. 

The accumulation of astaxanthin has been correlated 
with increasing fatty acid biosynthesis in H. pluvialis cells 
under stress conditions (Zhekisheva et al, 2002). In haema- 
tocysts, the fatty acid content can be as high as 30% of DW 
and comprised of >80% unsaturated fatty acids, making 



this organism a promising strain for biodiesel production 
(Zhekisheva et al., 2002; Miao et al, 2008; Damiani et al, 
2010; Peled et al, 2011; Lei et al, 2012). Under stress con- 
ditions, a large portion of de novo synthesized fatty acids 
are assembled into triacylglycerols (TAGs), which are co- 
sequestered along with astaxanthin esters into cytosolic 
lipid bodies through an unknown mechanism. Thus, unrav- 
elling the metabolic regulation of TAG biosynthesis may 
enable further elucidation of astaxanthin synthesis and 
deposition in H. pluvialis, which may in turn guide rational 
physiological or genetic manipulation to improve astax- 
anthin and lipid production in this organism. In the cur- 
rent study, a lipidomics approach was employed to analyse 
TAG molecular species and membrane glycerolipids quan- 
titatively in H. pluvialis. Lipidomics analysis will not only 
provide a broad context for understanding TAG biosyn- 
thesis in living cells, but will also help researchers under- 
stand membrane turnover and remodelling under high light 
acclimation. 

RNA sequencing (RNA-seq) is a newly developed, deep- 
sequencing technology for analysing entire transcriptomes. 
High-throughput DNA sequencing can produce a library 
of short cDNA reads that can then be aligned to a refer- 
ence genome or transcriptome, or can be assembled de novo 
without knowing the genome sequence of a given organism. 
This technology can provide a genome-scale transcription 
map that consists of the transcriptome structure and/or the 
level of transcriptional expression of individual genes (Wang 
et al, 2009b). RNA-seq is particularly attractive for unse- 
quenced organisms such as H. pluvialis because it does not 
rely on existing knowledge of the genome sequence for quan- 
tifying the transcriptome (Guarnieri et al, 2011; Rismani- 
Yazdi et al, 2011, 2012; Radakovits et al, 2012). 

In the current study, congruent lipidomic and transcrip- 
tomic analyses of H. pluvialis cells in two life cycle forms — 
macrozooids and haematocysts — were conducted to reveal its 
genomic and metabolic features. Special emphasis was placed 
on deciphering the regulation of astaxanthin biosynthesis 
and lipid metabolism during the encystment process under 
conditions of high irradiance. 

Materials and methods 

Cell growth conditions 

Haematococcus pluvialis (UTEX #2505) cells were cultured in glass 
columns (4cm diameter) containing MES-Volvox (MV) medium 
at 25 °C under continuous illumination (20 umol photon nr 2 s~'). 
Cultures were aerated with sterilized air containing C0 2 at 398 
ppm and at a feed velocity of 100 ml mhr 1 . Initial cell density was 
0.27 x 10 5 cells mr 1 . For high light level treatments, cells in the late 
exponential stage (2.23 x 10 5 cells ml -1 ) were exposed to continu- 
ous levels of high irradiance at 400 umol photon nr 2 s -1 for 48 h. 
Cell number was measured using a haemocytometer under light 
microscopy. Total chlorophyll contents were calculated using the 
method described by Porra (2002), and total carotenoid contents 
were measured using the calculation described by Lichtenthaler 
and Welburn (1983). For lipidomic analysis, cells were cultured 
under 400 umol photon mr 2 s for 4 d, and cells from days 0, 1,2, 
and 4 were harvested and lyophilized prior to lipid extraction and 
analysis. 
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Measurement of ROS production from H. pluvialis 

2', 7'-Dichlorodihydrofluorescein diacetate (H 2 DCFDA; Sigma, St 
Louis, MO, USA) was used to measure ROS according to the method 
described by Han et ol. (2012) with some modifications. A 5 ml aliquot 
of 1.35 x 10 5 cells ml -1 H. pluvialis was harvested by centrifugation at 
1000 g for 5min. Cells were resuspended in 5 ml of extraction buffer 
(50 mM TRIS-HC1, pH 7.5) supplemented with 5 uM H 2 DCFDA. 
The suspension was incubated at room temperature in darkness for 
lOmin. Algal cells were centrifuged at 1000 g for 5min at 4 °C and 
washed twice with 5 ml of extraction buffer. To extract the intracel- 
lular DCF, 200 ul of extraction buffer was added to the cell pellet 
and cells were lysed by 1 min of homogenization. An additional 4.8 ml 
of extraction buffer was added and centrifuged at 2000 g for 8 min 
to collect the supernatant. The supernatant liquid was diluted 1/100 
and the fluorescence was measured using a fluorescence spectropho- 
tometer LS55 (PerkinElmer, USA). The excitation and emission wave- 
lengths were 488 nm and 525 nm, respectively. Fluorescence values 
from untreated cells were calculated to correct the fluorescence values 
obtained from cells treated with H 2 DCFDA. Corrected fluorescence 
values were then divided by cell number. Macrozooid and haemato- 
cyst cells were harvested at 0, 1, and 2 d under high light irradiance. 



Pigment measurement 

For pigment extraction, 2 ml cell cultures were harvested by centrifu- 
gation (10 000 g, 2min) and supernatants were discarded. Cells were 
lysed by homogenization in 0.2ml of 90% (w/w) acetone for 1 min . 
Samples were centrifuged at 1 0 000 g for 2 min to obtain the superna- 
tants, which were then filtered through a 0.2 urn nylon filter. Filtrates 
were analysed on a Shimadzu Prominence high-performance liquid 
chromatography (HPLC) model LC-20AD equipped with a Waters 
Spherisorb S5 ODS2 cartridge column (4.6 x 250mm). HPLC analy- 
sis was conducted according to the method described by Park et di 
(2010). Pigments were identified by retention time and absorption 
spectra with reference to pigment standards (DHI 14C Centralen; 
Denmark). Concentrations of individual pigments were determined 
by HPLC using chlorophyll and carotenoid standards. 



RNA preparation 

Haematococcus pluvialis cells were frozen with liquid nitrogen and 
ground using a mortar and pestle. Total RNA was extracted using the 
Plant RNeasy Mini Kit (Qiagen, Valencia, CA, USA). Polyadenylated 
mRNA was isolated from 50 ug of total RNA samples using the 
FastTrack™ MAG Micro mRNA Isolation Kit (Life Technologies, 
Gaithersburg, MD, USA) according to the manufacturer's protocol. 
The isolated mRNA was quantified using an Agilent RNA 6000 
NANO Chip (Agilent, Wilmington, DE, USA). Subsequently, 200 ng 
of mRNA was used for constructing the cDNA library. 



cDNA library construction 

cDNA rapid library preparation began with fragmentation of 
mRNA using ZnCl 2 and heat. Cleaved RNA fragments were reverse 
transcribed into first-strand cDNA using reverse transcriptase and 
random hexamer primers (Roche, Switzerland). RNA templates 
were removed before a replacement strand was synthesized to gener- 
ate double-stranded cDNA. The double-stranded cDNA ends were 
polished (blunted) and short adaptors were ligated onto both ends. 
The adaptors provided priming sequences for both amplification and 
sequencing of the sample library fragments as well as the sequencing 
key, a short sequence of four nucleotides used by the system soft- 
ware for base calling. Purification with size selection was conducted 
using AMpure beads (Beckman Coulter, Brea, CA, USA). Lastly, 
the quality of the cDNA fragment library was assessed using a 2100 
BioAnalyzer (Agilent) and the library was quantified to determine 
the optimal amount for use as template for emulsion-based clonal 
amplification. 



Emulsion PCR 

Single 'effective' copies of template species from the cDNA library 
to be sequenced were hybridized to DNA Capture Beads. The 
immobilized library was then resuspended in amplification solution 
and the mixture was emulsified and amplified by PCR. After ampli- 
fication, the DNA-carrying beads were recovered from the emulsion 
and enriched. The second strands of the amplification products were 
melted away as part of the enrichment process, leaving the ampli- 
fied single-stranded DNA library bound to the beads. The sequenc- 
ing primer was then annealed to the immobilized, amplified DNA 
templates. 

454 FLX Titanium sequencing run 

Following amplification, the DNA-carrying beads were placed in 
the wells of a PicoTiterPlate (PTP) device, with each well holding 
single-stranded DNA-containing beads. The loaded PTP was then 
inserted into the Genome Sequencer FLX System and exposed to 
sequencing reagents. Information from all wells of the PTP was cap- 
tured simultaneously by a camera and was processed in real time by 
the onboard computer. 

Quantitative real-time PCR 

Total RNA was converted to cDNA using 5X reverse transcription 
master pre-mix (Elpis, Korea). cDNA was amplified using SYBR 
pre-mix (Takara, Japan) and a Thermal Cycler Dice Real Time 
System TP 8200 (Takara, Japan). Reverse transcription quantitative 
real-time PCR (RT-qPCR) conditions were as follows: 95 °C for 10 
s; 95 °C for 5 s, 55-65 °C for 30 s (40 cycles); 95 °C for 15 s, 60 °C for 
30 s, 95 °C for 15 s; 72 °C for 5 min; and 15 °C for 5 min. To calculate 
relative changes in gene expression, results were analysed using the 
2-AA CT method described by Livak and Schmittgen (2001). 

All RT-qPCR experiments described were performed at least twice 
with three technical repeats (n=3) for independently prepared biological 
samples («=2). Results were represented as the mean ±SD from six meas- 
urements. The actin gene was used as an internal control, as it has been 
used in other study (Eom et ol. , 2006). The ratio between the expression 
values of each gene and the actin gene (Actl) in Haematococcus cells 
under normal light conditions (control) was set to 1 and the expression 
ratios of other conditions were given relative to the control. 

Lipid extraction from H. pluvialis 

For lipidomics analysis, a lOmg biomass was used for lipid extrac- 
tion. Total lipids were extracted by grinding the cell pellet using an 
ice-cold mortar and pestle in 6ml of chloroform:mefhanol (2:1, v/v). 
Lipid extracts were vigorously agitated at room temperature for 1 h. 
Extracts were then mixed with 1.5ml of potassium chloride (0.7%, 
w/v). Following centrifugation at 1000 g for 5 min, extracts were 
separated into two phases and total lipids in the lower organic phase 
were transferred into another vial using a glass pipette. The organic 
solvent was evaporated under a stream of nitrogen. 

Lipidomic analysis 

Total lipids were recovered in 1ml of chloroform:methanol (1:1). 
Lipidomic analysis was performed on a 6460 triple quadruple elec- 
trospray ionization mass spectrometer (ESI/MS) coupled with a 1260 
high-performance liquid chromatograph (Agilent). Single-stage mass 
spectrometry (MS) and MS/MS analyses in the positive or negative 
ion mode were conducted to identify lipid molecular species. Based 
on the fragmentation of lipid molecules, multiple reaction monitor- 
ing was developed and used for quantification. A 10 aliquot of lipid 
extract was mixed with internal standards and chloroform:methanol 
(1:1) to a final volume of 50 ul. For positive ion mode analysis, 
the internal standards mixture provided 2 uM digalactosyldiacyl- 
glycerol (DGDG; 18:0/18:0), 5 uM monogalactosyldiacylglycerol 
(MGDG; 18:0/18:0), 0.25 uM phosphatidylcholine (PC; 18:1/18:1), 
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0.32 uM phosphatidylethanolamine (PE; 17:0/14:0), and 1 uM TAG 
(17:0/17:0/17:0). For negative ion mode analysis, the mixture con- 
tained internal standards providing 1.243 uM phosphatidylglycerol 
(PG; 17:0/20:4) and 1.15 uM phosphatidylinositol (PI; 17:0/20:4). 
Prior to infusion into the ESI/MS, 1 ul of lipid mixture was sepa- 
rated at 40 °C on a ZOBAX SBC18 column (1.8 um, 2.1 x 150mm, 
Agilent) and an Extend C18 column (1.8 um, 2.1 x 150mm, Agilent) 
for positive and negative mode MS analyses, respectively. For positive 
mode, the mobile phases were methanol:acetonitrile:H 2 0 (19:19:2; 
A) and isoproponal (B) containing 0.1% formic acid and lOmM 
ammonium acetate. The liquid chromatography (LC) gradients 
were as follows: Omin, 90% A and 10% B; 5min, 90% A and 10% B; 
25min, 60% A and 40% B; 60min, 45% A and 55% B; 66min, 45% 
A and 55% B; and 68min, 90% A and 10% B. For negative mode, the 
mobile phases were 85% methanol (A) and isopropanol (B) contain- 
ing 0.025% NH 4 OH. The LC gradients were as follows: Omin, 95% 
A and 5% B; 15min, 85% A and 15% B; 22min, 45% A and 55% B; 
42min, 45% A and 55% B; and 44min, 95% A and 5% B. The flow 
rate was 0.2 ml min -1 . Lipid standards for the membrane glycerolipids 
were purchased from Avanti Polar Lipids (Alabaster, AL, USA) and 
TAG standards were purchased from Sigma-Aldrich. For lipidom- 
ics analysis, three biological replicates (w=3) were included. For each 
biological replicate, lipids were extracted in duplicate for independent 
LC/MS analysis. Lipid data were presented as the mean ±SD from 
six measurements for each time point. 



Results and Discussion 

Morphological and biochemical changes during 
encystment under high irradiance 

The life cycle of H. pluvialis consists of four distinct types 
of cells: microzooids, motile macrozooids, non-motile 



palmelloids, and haematocysts (Elliot, 1934). Under favour- 
able conditions, H. pluvialis cultures were dominated by mac- 
rozooids (Fig. 1A); under conditions of high light intensity 
(400 umol photon nr 2 s" 1 ), 2 d cultures of H, pluvialis were 
transformed into haematocysts while accumulating large 
amounts of carotenoids (Fig. IB; Supplementary Table 
S2 available at JXB online). The total carotenoid content 
increased from 1.5 |ig ml" 1 to 6.22 |ig mL 1 after 2 d of expo- 
sure to high irradiance conditions, while total chlorophyll 
contents were reduced from 2.32 itg mL 1 to 1.60 |ig mT 1 
(Supplementary Table S2 available at JXB online). In green 
macroozoid cells, chlorophyll a and b and lutein were the 
most abundant pigments, followed by (3-carotene. In contrast, 
astaxanthin acyl-esters and (3-carotene were the prevailing 
pigments in haematocysts, followed by lutein and chlorophyll 
a and b. As indicated by fluorescence emitted after Nile red 
staining, cells exposed to high light formed a large number of 
carotenoid-rich lipid bodies in the cytosol, a striking oleagi- 
nous characteristic of haematocysts (Fig. 1C). 

To investigate cellular alterations in glycerolipid compo- 
sition and content during encystment, glycerolipidomes of 
Haematococcus cells exposed to high light were analysed 
quantitatively. The total content of glycerolipids increased 
from 75.99 umol g _1 DW to 96.80 umol g" 1 DW following 
exposure of H. pluvialis cells to high light for 4 d (Fig. 2A). 
Glycolipids including MGDG, DGDG, and sulphoquinovo- 
syl diacylglycerol (SQDG), the main constituents of the chlo- 
roplast membranes and photosynthetic complexes, were the 
major lipid classes in macrozooid cells (day 0) and accounted 
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for up to 67.8% of total glycerolipids. Among the extraplas- 
tidic membrane glycerolipids including diacylglyceryltrimeth- 
ylhomoserine (DGTS), PC, PI, and PE, DGTS was the most 
abundant (~5.6% of total glycerolipids). The presence of two 
structural analogues, DGTS and PC, distinguished H. pluvia- 
lis from Chlamydomonas reinhardtii, which does not contain 
PC in its lipidome (Yoon et al, 2012). After 4 d under high 
light, TAG in macrozooid cells increased from 10.26 umol 
g" 1 DW (13.5% of total glycerolipids) to 65.91 umol g 1 DW 
(67.8% of total glycerolipids). A concomitant and remarkable 
decrease of MGDG, DGDG, and SQDG and slight reduc- 
tions of PE, DGTS, and PG were also observed, while PC 
and PI remained stable (Fig. 2B). For example, MGDG and 
DGDG accounted for 21.9% (corresponding to 16.64 umol 



g" 1 DW) and 17.6% (13.35 umol g" 1 DW) of total glycerolipids 
in green macrozooid cells, and decreased to 6.3% (6.11 umol 
g" 1 DW) and 4.7% (4.52 umol g" 1 DW) of total glycerolip- 
ids, respectively, in haematocyst cells after 4 d of high light 
induction. SQDG decreased from 21.53 umol g" 1 DW (28.3% 
of total glycerolipids) to 11.73 umol g" 1 DW (12.1% of total 
glycerolipids) after 4 d under high irradiation (Fig. 2A). 

As shown in Fig. 3, high irradiance induction exerted a 
profound impact on the cellular content of individual TAG 
species and molecular species of glycolipids. The major TAG 
species (52:2, 52:4, 52:5, 52:8, and 50:1) accounted for >50% 
of total TAG, among which TAG 52:2 was the most predomi- 
nant, accounting for 24.7% of total TAG (Fig. 3A). Eight spe- 
cies (50:2, 50:3, 50:9, 52:1, 52:6, 52:7, 52:9, and 54:10) were 
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designated as the minor TAG species, which accounted col- 
lectively for 20-40% of total TAG. The other 11 TAG species, 
accounting for <5% of total TAG, were defined as trace TAG 
molecular species (Supplementary Fig. SI available at JXB 
online). All TAG species increased in response to high irradi- 
ance conditions. 

Four species of MGDG were detected, namely 34:4, 34:5, 
34:6, and 34:7. MGDG 34:4, 34:5, and 34:6 decreased dra- 
matically after 1 d under high irradiance and then decreased 
more moderately from day 2 to day 4, whereas 34:7 remained 
at a steady level over 4 d (Fig. 3B). In contrast to MGDG, the 
molecular species of the bilayer glycerolipid DGDG exhibited 
complex responses to high irradiance induction. For example, 
among four DGDG species, DGDG 34:2 and 34:4 progres- 
sively decreased over 4 d under high irradiance conditions, 34:3 
slightly increased on day 1 and then declined to a level lower 
than the basal level (day 0), and DGDG 34:5 was not affected 
during the first 2 d and then decreased by 3-fold in the follow- 
ing 2 d (Fig. 3C). Three major SQDG species (32:0, 34:2, and 
34:3) decreased substantially when induced by high irradiance 
conditions, while the minor species SQDG 34:4 decreased only 
slightly (Fig. 3D). The profiles of PG, PC, DGTS, PI, and PE 
from Haematococcus cells exposed to high irradiation condi- 
tions are presented in Supplementary Fig. S2 available at JXB 
online. 

Sequencing and de novo assembly of the H. pluvialis 
transcriptome 

The cDNA libraries for macrozooids and haematocysts were inde- 
pendently prepared twice to rrrinimize any biological variations 



Four pyrosequencing runs using a Roche GS-FLX 454 yielded a 
total of 2 141 048 raw reads, equivalent to -848 Mb of data for the 
H. pluvialis transcriptomes The average length of the raw sequenc- 
ing reads was ~396bp. The pre-assembly process of trirrrming and 
cleaning for quality control yielded a total of 2 139 71 1 high qual- 
ity reads, which were further assembled into 1 453 995 contiguous 
sequences (contigs). A total of 17 364 reads remained as singletons 
(reads not incorporated into the assembly). The average length 
of contigs was 589 bp, and overlapping contig reads were further 
assembled into 34 947 isotigs. The average length of the isotigs was 
1555 bases with an N50 of 1737 bases, meaning that 50% of the 
assembled bases were incorporated into isotigs >1737 bases A sta- 
tus summary of the sequencing and assembly results is presented 
in Table 1 and the size distributions for raw reads, assembled con- 
tigs, isotigs, and sequence coverage are presented in Supplementary 
Fig. S3 available at JXB online. Isotigs sharing common contigs 
were clustered into individual isogroups, with each presumably 
representing one gene locus capable of producing multiple, alterna- 
tively spliced transcripts. Lastly, the de novo transcriptome assem- 
bly identified a total of 18 482 isogroups, possibly representing the 
total number of genes in the H, pluvialis transcriptome under the 
given culture conditions. The 6937 isogroups contained multiple 
isotigs (equivalent to an average of 1.89 transcripts per gene), while 
1 1 545 isogroups had only one isotig or contig each. Raw data files 
are available at the National Center for Biotechnology Information 
(NCBI) Sequence Read Archive (accession no. SRR1040551). 

Functional genome annotation 

All unique sequences (total of 267 contigs and 34 680 isotigs) 
were queried against the sequences available in the NCBI 
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Table 1 . Summary of de novo assembly in H. pluvialis 



Sequencing 


No. of sequences 


Bases 


Raw sequencing reads 


2 141 048 


848 713 419 


Average read length 




396 


Pre-assembly 






Trashed 


1337 


3 795 275 


Reads used in assembly 


2 139 711 


844 918 144 


Average read length 




395 


Assembly 






Contigs 






Reads assembled as contigs 


1 453 995 


558 385 452 


Number of contigs 


62 636 


36 715 613 


Average length of contigs 




589 


Isotigs 






Number of isotigs 


34 947 


54 330 087 


Average length of isotigs 




1555 


Isotig N50 




1737 


Number of isogroups 


18 482 




Singletons 


17 364 





non-redundant protein database using the BLASTx algo- 
rithm for annotation. Approximately 57.9% of the BLAST 
queries (20 234 unique sequences) produced significant blast 
matches with expected values (E-values) below a threshold 
of 10~ 6 . The top-hit species distribution of BLAST matches 
revealed distinct characteristics of the microalgal transcripts. 
For example, 70% of the total BLAST matches were clustered 
to related microalgal species, predominantly Volvox carteri 
f. nagariensis (39%) and C. reinhardtii (31%) (Fig. 4A). 

The best hits (sequences producing the highest scores/low- 
est E-values) selected from the BLAST queries were further 
annotated with gene ontology (GO) terms through the GO 
site (http://www.geneontology.org/). The GO database rep- 
resents ontological descriptions for gene products according 
to three main categories: biological processes, cellular com- 
ponents, and molecular functions. A total of 61 GO terms 
were assigned to 22 431 cDNAs, including isotigs and con- 
tigs. Among all GO terms extracted, 7959 (23%) belonged to 
the molecular functions class and 7291 (21%) and 7181 (20%) 
were grouped into the biological processes and cellular com- 
ponents classes, respectively (Fig. 4B). 

Of the assembled isotigs and contigs, 36% were not sub- 
ject to standardized GO classification because they did not 
match any sequence in the GO database. Such a large quan- 
tity of 'no hit' sequences may have been attributable to a 
lack of transcriptome data or gene prediction data of simi- 
lar species in the GO database. Of the remaining assembled 
isotigs and contigs, a significant number of sequences were 
associated with biological processes referred to as metabolic 
process (26% of GO-annotated sequences searched), fol- 
lowed by those that were unclassified (16%), those referred 
to as cellular process (12%), and those associated with bio- 
logical regulation (11%). As for those sequences associated 
with molecular functioning, genes related to binding (39%) 
and catalytic activity (35%) were the most abundant, fol- 
lowed by those that were unclassified (12%). In the cellular 



components category, ~39% of the transcripts belonged to 
cell part, 12% were unclassified, 12% belonged to organelle, 
and 9% to membrane. A detailed analysis of GO terms of the 
three function classes is presented in Supplementary Fig. S4 
available at JXB online. 

Kyoto Encyclopedia of Genes and Genomes (KEGG) 
pathway analysis is an alternative approach to categorizing 
gene functions with an emphasis on biochemical pathways. To 
reconstruct the metabolic pathways in H. pluvialis, all assem- 
bled unique sequences were queried against the KEGG data- 
base and 22.4% of queries were assigned to KEGG orthology 
(KO) identifiers (ID) with E-values below the threshold value 
of 10~ 6 . These KO IDs were mapped against all pathway maps 
through the basic pathway mapping tool of the KEGG map- 
per (http://www.genome.jp/kegg/mapper.html). As shown in 
Fig. 5, H. pluvialis expressed transcripts for enzymes func- 
tioning in most of the major metabolic pathways, including 
metabolism of N- or O-glycans, lipids, carbohydrates, energy, 
amino acids, nucleotides, cofactors and vitamins, and isopre- 
noids. However, transcriptional expression of most of the 
genes related to certain pathways — such as the biosynthetic 
pathways for glycosphingolipids, phenylpropanoids, and 
alkaloids, and xenobiotic biodegradation and metabolism — 
were not detected. 

Comparison of global transcriptomes of macrozooids 
and haematocysts 

Normalized abundances of individual transcripts (assem- 
bled sequences including isotigs and contigs) in green mac- 
rozooids and haematocysts are presented in a scatter plot in 
Fig. 6. Over 2-fold differential expression was detected in as 
many as 38.1% of the assembled sequences (13 313 out of 
34 947) between the two cell types. Surprisingly, ~30.6% of 
the transcripts (10 705/34 947) were expressed more abun- 
dantly in haematocysts than in macrozooids. In contrast, 
only 7.5% (2608/34 947) of the transcripts were expressed in 
higher amounts in macrozooids than in haematocysts. As two 
independently prepared cDNA library sets were available, the 
same analysis was performed separately for each set, confirm- 
ing that haematocysts up-regulated many more genes than 
green macrozooid cells. Multiple lines of evidence indicated 
specificity regarding the high irradiance-induced up-regula- 
tion of many genes of great diversity (as much as ~30% of all 
unique sequences analysed). First, the alteration pattern of 
many sets of genes revealed good agreement with the known 
physiology of Haematococcus cells under high light condi- 
tions. For example, the coordinated and apparent increase in 
transcript abundance was noticeable in many genes related 
to protection against photooxidative stress. Conversely, the 
expression of a number of highly active genes in macrozoo- 
ids, especially genes encoding proteins involved in light-har- 
vesting activity, declined under high irradiance conditions. 
Secondly, the extent of up-regulation varied among different 
pathways/genes. For example, the de novo fatty acid biosyn- 
thesis pathway was up-regulated to a greater extent than the 
TAG assembly pathway under high irradiance conditions. 
Thirdly, the transcriptome data suggested that >61% of 
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Fig. 4. Top-hit species distribution of BLAST matches and gene ontology (GO) analysis of H. pluvialis. (A) Total BLAST matches were clustered to related 
green microalgae species, Volvox carted f. nagariensis (39%) and Chlamydomonas reinhardtii (31 %). (B) Isotigs and contigs of H. pluvialis were categorized 
in three classes: molecular function (23%), biological process (21 %), and cellular component (20%). (This figure is available in colour at JXB online.) 
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Fig. 5. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of H. pluvialis transcripts. Isotigs and contigs of expected values (E-values) 
below the threshold of 1 0" 6 were analysed. Haematococcus pluvialis expressed transcripts for enzymes functioning in most of the major metabolic 
pathways, including N- or O-glycan, lipid, carbohydrate, energy, amino acid, nucleotide, cofactor and vitamin, and isoprenoid metabolism. 



genes were not responsive to high irradiance, such as those housekeeping genes in the next-generation sequencing (NGS) 
involved in the cell cycle, carbon fixation, and chlorophyll remained steady in their expression under high irradiance 
metabolism, etc. Lastly, most of the multiple representative (Supplementary Table S3 available at JXB online). In fact, 
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Fig. 6. Scatter plot of normalized abundances of each cell type [green 
(macrozooid) and red (haematocyst)]. Differential expression >2-fold was 
detected in as many as 38.1 % of assemblages between macrozooid and 
haematocyst cells. Approximately 30.6% of the transcripts were expressed 
more abundantly in haematocyst cells than in macrozooid cells. (This figure 
is available in colour at JXB online.) 



the expression of >30 different kinds of genes which have 
been used as housekeeping gene controls in various studies 
was examined between macrozooids and haematocysts in 
the NGS data, and the slope value (0.9372) of the correla- 
tion graph demonstrated that the expression levels of each 
gene between the two kinds of cells was quite similar overall 
(Supplementary Fig. S5 available at JXB online). 

The assembled sequences that showed differential expres- 
sion (>2-fold) between haematocysts and macrozooids 
were mapped on the KEGG reference pathway diagram 
(Supplementary Fig. S6 available at JXB online). The KEGG 
analysis revealed that those sequences up-regulated in haema- 
tocysts included various genes encoding proteins involved in a 
wide range of metabolic processes in addition to the genes for 
the synthesis of astaxanthin and oil, the two major products 
of haematocysts. For example, many genes coding for vari- 
ous metabolic enzymes, particulary for the biosynthesis of 
secondary metabolites, microbial metabolism in diverse envi- 
ronments, and nucleic acid and amino acid metabolism were 
more highly expressed in Haematococcus cells under high 
irradiance conditions. Additionally, genes up-regulated in 
haematocysts were enriched with those known to be involved 
in processes related to gene expression, such as spliceosome 
formation, RNA transport, ribosome biogenesis, and protein 
processing. Such an increase in abundance of those transcripts 
related to the information flow may aid the active translation 
of the transcripts of the diversity under high light condi- 
tions. The present NGS data analysis of haematocyst tran- 
scripts also indicated that most of the genes involved in basic 
metabolism such as glycolysis, the tricarboxylic acid (TCA) 
cycle, and oxidative phosphorylation were expressed at levels 
comparable with those of macrozooids and that few genes 



in those pathways were even up-regulated. Collectively, this 
comparative transcriptomic analysis suggested that haemato- 
cysts, previously referred to as resting cells, were probably not 
quiescent at all, but instead were probably more metabolically 
active than macrozooids. 

Analysis of the transcriptome involved in astaxanthin 
biosynthesis 

Multiple lines of transcriptome evidence indicated that syn- 
thesis of the five-carbon isoprene units in H. pluvialis was 
most probably mediated by the MEP pathway rather than 
the MVA pathway (Supplementary Fig. S7 available at JXB 
online). For example, transcripts for all eight genes involved 
in the MEP pathway were identified in Haematococcus cells. 
Four of the MEP pathway genes exhibited higher transcript 
levels in haematocysts — as much as ~1.8- to 2-fold of that 
in macrozooids — which correlated with the accumulation of 
astaxanthin. The four up-regulated genes were those encoding 

1- deoxy-D-xylulose-5-phosphate synthase (DXS), 1-deoxy- 
D-xylulose 5-phosphate reductoisomerase (DXR), 1-hydroxy- 

2- methyl-2-(E)-butenyl-4-diphosphate synthase (HDS), 
and 4-hydroxy-3-methylbut-2-enyl diphosphate reductase 
(HDR), which showed the highest fold change. In contrast, 
only two types of MVA pathway genes, namely one encod- 
ing acetyl-CoA C-acetyltransferase (ACAT) and two encod- 
ing hydroxymethylglutaryl-CoA (HMGS), were identified in 
the transcriptome of either haematocysts or macrozooids. Of 
these putative MVA pathway genes, only one of the HMGS 
genes was up-regulated at the transcript level, whereas the 
others did not change under high irradiance conditions. 

DMAPP synthesis is thought to be catalysed by isopente- 
nyl diphosphate isomerase (IPPI) (Sun et al, 1998). However, 
the transcriptomic analysis indicated that neither of the two 
putative IPPI transcripts was up-regulated during the cel- 
lular accumulation of astaxanthin (Fig. 7). The results of 
previous studies have suggested that HDR was an additional 
enzyme capable of catalysing interconversion between IPP 
and DMAPP (Hoeffler et al, 2002; Rohdich et al, 2002). 
One of the two HDR-encoding transcripts was moderately 
up-regulated (>2-fold) under high irradiance (Supplementary 
Fig. S7 available at JXB online), indicating that this enzyme 
was more likely to be responsible for supplying the isoprene 
units for astaxanthin synthesis under high light conditions in 
this organism. 

The astaxanthin biosynthetic pathway from IPP and DMAPP 
is presented in Fig. 7. All of the genes essential for astaxanthin 
synthesis from DMAPP were up-regulated in concert under 
high irradiance with the exception of IPPI. As the first com- 
mitted step in the pathway, phytoene biosynthesis has long been 
considered a bottleneck in carotenoid biosynthesis (Lorenz and 
Cysewski, 2000; Cardozo et al. , 2007). Previously, the expression 
of the phytoene synthase (PSY) gene in Haematococcus cells 
was shown to be up-regulated by high light levels (Steinbrenner 
and Linden, 2003; Vidhyavathi et al, 2008). Of the four phy- 
toene synthase cDNA sequences, including partial sequences 
reported for H. pluvialis in GenBank, the present analysis con- 
firmed that only one was up-regulated in haematocysts. The 
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Fig. 7. Differential expression of genes involved in the astaxanthin biosynthesis pathway. Coloured circles indicate the normalized read number ratio of 
red (haematocyst) per green (macrozooid) cells and each circle represents an independent isogroup (gene). The star designates the primer-designed 
gene for RT-qPCR. The bar graph indicates the expression pattern of each star-marked gene in the RT-qPCR analysis. Values in the bar graphs are 
the mean ±SD of six different measurements. IPPI, isopentenyl diphosphate isomerase; FPPS, farnesyl diphosphate synthase; GGPS, geranylgeranyl 
diphosphate synthase; PSY, phytoene synthase; PDS, phytoene desaturase; ZDS, zeta-carotene desaturase; LBC, lycopene beta-cyclase; CHY, 
carotenoid hydroxylase; BKT, beta-carotene ketolase. 



gene encoding phytoene desaturase (PDS) that catalyses the 
first step in lycopene biosynthesis was also up-regulated at the 
mRNA level under high irradiance (Steinbrenner and Linden, 
2003; Steinbrenner and Sandmann, 2006). The transcriptome 
analysis demonstrated that two PDS genes were up-regulated by 
~1.5- to 1.9-fold in haematocysts compared with the macrozoo- 
ids. The plastid terminal oxidase (PTOX) has been suggested to 
serve as a cofactor for desaturation of phytoene and (3-carotene 
(Li et al, 2008). Previously, two PTOX genes (PTOX1 and 2) 
have been reported in H. pluvialis (Wang et al. , 2009a; Houille- 
Vernes et al, 2011), and our RNA-seq analysis identified their 
transcripts. As shown in Table 2, transcripts of the PTOX1 
and 2 genes increased >2-fold under high light level condi- 
tions. Consistent with previous observations (Steinbrenner and 
Linden, 2003), the transcripts of the Haematococcus lycopene 
|3-cyclase gene were >2-fold more abundant in haematocysts 
than in macrozooids. It is worth noting that expression of the 
e-cyclase and zeaxanthin epoxidase genes was either undetect- 
able or slightly reduced, suggesting that lycopene is largely 
shunted into (3-carotene synthesis as part of astaxanthin syn- 
thesis. Subsequent oxygenation reactions converting (3-carotene 
to astaxanthin are considered rate-limiting steps catalysed by 
(3-carotene ketolase (BKT) and |3-carotene hydroxylase (CHY) 
(Steinbrenner and Linden, 2001, 2003; Vidhyavathi et al. , 2008). 
The present NGS data showed that these two genes were also 
up-regulated in haematocysts. 

Collectively, almost all of the genes involved in astaxan- 
thin biosynthesis from IPP were found to be simultaneously 



up-regulated in transcriptional expression under high irradi- 
ance (Fig. 7). The concerted reactions of the multiple enzy- 
matic steps in carotenogenesis have already been demonstrated 
as critical for the enhancement of astaxanthin biosynthesis. 
For instance, the introduction of yeast astaxanthin biosynthe- 
sis genes in Saccharomyces cerevisiae successfully enhanced 
the accumulation of astaxanthin and its intermediates in 
transgenic cells (Ukibe et al, 2009). Up-regulation of the 
carotenogenic genes responsible for astaxanthin biosynthesis 
from GGPP also led to excessive astaxanthin accumulation 
in a mutant strain of Phaffia rhodozyma (Miao et al, 2011). 
Therefore, the greater transcript abundance of the multiple 
genes involved in astaxanthin biosynthesis may contribute 
collectively to the observed increase in astaxanthin synthe- 
sis under high irradiance. In fact, the additional RT-qPCR 
analysis revealed that up-regulation of all genes involved in 
astaxanthin biosynthesis occurred during the first 24 h under 
high irradiance. The most dramatic increases in transcript 
abundance (up to 20- to 30-fold) were observed in genes such 
as PSY, PDS, and CHY. The high level of responsiveness of 
these genes was considered indicative of their prominent regu- 
latory roles in augmenting astaxanthin biosynthesis. 

Transcriptome analysis related to lipid metabolism in 
H. pluvialis 

Under stress conditions, many green microalgae can increase 
the number and size of their lipid bodies, which are discrete 
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Table 2. Isotigs ofROS scavenging-related genes in the H. pluvialis transcriptome 



Isogroup 


Isotig number 


Description 


Accession 


Read ratio (red/green) 


isogroup04416 


isotigl 81 19-18120 


Catalase [Haematococcus pluvialis) 


ABK41 476.1 


1.40 


isogroup00630 


isotig0731 7-07322 


Mn superoxide dismutase 4 (Chlamydomonas 
reinhardtii) 


ACZ37254.1 


1.07 


isogroup05851 


isotig20989-20990 


Mn superoxide dismutase (Haematococcus 
pluvialis) 


AAW69292.1 


1.08 


isogroup05972 


isotig21231-21232 


Mitochondrial Mn superoxide dismutase 
(Haematococcus pluvialis) 


ABF82431.1 


1.67 


isogroup11662 


isotig28274 


Chloroplast Fe superoxide dismutase 1 precursor 
(Chlamydomonas reinhardtii) 


XP_001 690591.1 


1.15 


isogroup04801 


isotigl 8889-1 8890 


L-Ascorbate peroxidase (Volvox carter! 
f. nagariensis) 


XP_002957861.1 


1.78 


isogroupl 1 649 


isotig28261 


L-Ascorbate peroxidase (Chlamydomonas 
reinhardtii) 


XP_001 695476.1 


2.37 


isogroup00937 


isotig08842-08845 


Glutathione reductase (Coccomyxa subellipsoidea 
C-169) 


EIE1 8555.1 


1.48 


isogroup04339 


isotigl 7965-1 7966 


Glutathione reductase (Chlamydomonas reinhardtii) 


XP_001 696579.1 


1.58 


isogroup05030 


isotigl 9347-1 9348 


Chloroplast glutathione reductase (Solanum 
lycopersicum) 


NP_001 234243.1 


3.82 


isogroup0401 0 


isotigl 7307-1 7308 


Plastid terminal oxidase (Haematococcus pluvialis) 


ABV72392.1 


2.11 


isogroup04343 


isotigl 7973-1 7974 


Plastid terminal oxidase (Haematococcus pluvialis) 


ABV72391.1 


2.07 


isogroup08668 


isotig25280 


Monodehydroascorbate reductase (Coccomyxa 
subellipsoidea C-169) 


EIE20462.1 


0.79 


isogroup21067 


isotig37679 


Dehydroascorbate reductase (Chlamydomonas 
reinhardtii) 


XP_001 698375.1 


10.84 



TAG- and carotenoid-containing organelles surrounded 
by a lipid monolayer (Thompson, 1996; Boussiba, 2000; 
Grunewald et al., 2001). During cyst formation under stress 
conditions, H. pluvialis accumulates large amounts of TAG 
in extraplastidic lipid bodies, accompanied by astaxanthin 
deposition. Results from previous studies of the unicellular 
microalga Dunaliella bardawil demonstrated that p-carotene 
accumulation was closely correlated with TAG synthesis and 
oil globule formation (Katz etal., 1995; Rabbani etal, 1998). 
Such sink activity of oil globules for secondary carotenoids 
is assumed to be critical for the massive accumulation of 
pigments in algal cells, as it may prevent end-product inhibi- 
tion of the carotenoid biosynthetic pathway (Rabbani et al. , 
1998). Zhekisheva et al. (2002) previously found that H. plu- 
vialis actively synthesized fatty acids de novo — mainly oleic, 
palmitic, and linoleic acids — and deposited them primarily in 
the form of TAGs and astaxanthin esters. Accumulation of 
astaxanthin in these algal cells was also positively correlated 
with that of newly synthesized oleic acid (Zhekisheva et al, 
2002). 

Figure 8 shows the microalgal fatty acid synthetic process 
in the chloroplast. In the first committed step of fatty acid 
synthesis, acetyl-CoA is converted to malonyl-CoA precur- 
sors by acetyl-CoA carboxylase (ACCase). This is followed 
by a series of successive condensation reactions, ultimately 
resulting in the production of an acyl-acyl carrier protein 
(ACP). The transfer of a malonyl group from malonyl-CoA 
to the ACP is catalysed by malonyl-CoA:ACP transacylase 
(MCAT), forming malonyl-ACP. An acyl-ACP chain is suc- 
cessively elongated by a subsequent series of condensation 



reactions catalysed by 3-ketoacyl-ACP-synthase (KAS), 
3-ketoacyl-ACP-reductase (KAR), 3-hydroxyacyl-ACP 
dehydratase (HD), and enoyl-AC reductase (ENR), yielding 
C16:0- and C18:0-ACP. In the de novo fatty acid biosynthesis 
process, a double bond can be introduced to the A9 position 
of C18:0-ACP by a stearoyl-ACP desaturase (SACPD) to 
form CI 8:1 -ACP. Termination of elongation can be medi- 
ated either by an acyl-ACP thioesterase (FAT) — resulting 
in free fatty acid release and export to the cytosol — or by 
direct transfer of the acyl group to glycerol-3-phosphate and/ 
or monacylglycerol-3-phosphate in the TAG biosynthetic 
pathway (Ohlrogge and Browse, 1995; Hu et al, 2008). In 
the transcriptome analysis related to fatty acid biosynthesis, 
all transcripts with the exception of ENR were identified, 
including single genes for MCAT, KAS I, KAS III, KAR, 
HD, SACPD, and FAT, two genes for KAS II and ACP, and 
six genes for ACCase. Among those genes, the levels of KAR 
and KAS I transcripts were increased to the greatest degree 
(2.75- and 3.5-fold, respectively) in response to conditions 
of high light. Expression levels of ACCase, MCAT, KAS II, 
III, HD, and SACPD were similar between haematocysts and 
macrozooids. However, RT-qPCR analysis revealed that the 
fatty acid biosynthetic genes, with the exception of HD, were 
up-regulated at least 2- to 15-fold in their expression during 
24 h of high irradiation conditions (Fig. 8). Collectively, the 
results of the gene expression analyses suggested that high 
irradiance stimulated fatty acid biosynthesis in either a tran- 
sient or progressive manner in H. pluvialis. 

Multiple pathways for TAG biosynthesis have been elu- 
cidated in eukaryotes. The acyl-CoA-dependent pathway 
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Fig. 8. Differential expression of genes involved in the fatty acid biosynthesis pathway. Coloured circles indicate the normalized read number ratio of red 
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desaturase. 



involves three sequential transfers of a fatty acyl moiety from 
acyl-CoA, the last of which is acylation of a diacylglycerol 
(DAG) with an acyl-CoA to form TAG. The first fatty acid 
transfer to position 1 of glycerol-3-phosphate (G3P) produces 
lyso-phosphatidic acid (LPA) by the catalysis of G3P acyl- 
transferase (GPAT). Another acyl moiety transfer to posi- 
tion 2 of LPA results in the formation of phosphatidic acid 
(PA), which is catalysed by LPA acyltransferase (LPAAT). 
Formation of DAG from dephosphorylation of PA is facili- 
tated by phosphatidic acid phosphatase (PAP), and the final 
transfer of a third acyl chain to position 3 of DAG is acom- 
plished by diacylglycerol acyltransferase (DGAT), leading 
to synthesis of a neutral TAG. In the present NGS analy- 
sis, a single gene was identified for GK, LPAAT, and PAP. 
Additionally, two and three genes were identified for GPAT 
and DGAT, respectively. While the transcript levels of GPAT 
were very similar between macrozooids and haematocysts, 
those for LPAAT, PAP, and DGAT were expressed at slightly 
higher levels (~1.5- to 1.9-fold) in the red cysts (Fig. 9). 



Additionally, TAG can be synthesized through an acyl-CoA- 
independent pathway mediated by a phospholipid:DAG 
acyltransferase (PDAT) in microalgae and other eukary- 
otes. PDAT transfers a fatty acyl group from a phospholipid 
(PL) or galactolipid to DAG, leading to production of TAG. 
Additionally, PDAT has been suggested to function as an 
acyl-CoA-independent diacylglycerol transacylase (DGTA) 
in C. reinhardtii, which produces TAG by using two DAG 
molecules as both acyl donors and acceptors (Yoon et ah, 
2012). In the present analysis, a single gene was identified for 
PDAT, which was slightly up-regulated (1.6-fold) under high 
irradiance. Early expression analysis (0-24 h) via RT-qPCR 
revealed that the expression of LPAAT was highly up-regu- 
lated (10-fold), and that of all other genes except DGAT were 
moderately up-regulated during 24 h under high irradiance 
(Fig. 9). 

Oil globule-associated proteins are known to be involved 
in various processes, such as globule formation, degrada- 
tion, stabilization, and in globule-globule interactions or 
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interactions between globules and other organelles (Katz 
et al, 1995; Frandsen et al, 2001; Murphy, 2001, 2004). Three 
genes encoding the major oil globule-associated protein, the 
structural component of the monolayer of oil globules (Peled 
et al. , 201 1), did not show any differences between the cells at 
the two different stages (data not shown). This result implied 
that the globule-associated protein of this organism might 
be regulated at the translational or post-translational level to 
sustain the biogenesis and enormous accumulation of the oil 
globules. 

Although the cellular content of the membrane glyc- 
erolipids substantially decreased or remained unchanged in 
response to elevated light intensity, the transcriptome analy- 
sis revealed that a number of genes involved in membrane 
glycerolipid biosynthesis and turnover were up-regulated at 
the gene expression level (Supplementary Table S4 available 
at JXB online). Lipidome and transcriptome results taken 
together allowed for speculation that the de novo biosynthesis 



of the membrane glycerolipids was probably enhanced in 
response to the increased irradiance and therefore provided 
newly synthesized membrane glycerolipid molecules with 
the proper structure and function to replace damaged and 
degraded lipid species in the cellular membrane or photosyn- 
thetic complexes. 

Transcriptome analysis of additional photoprotective 
mechanisms under high irradiance 

Over-reduction of the photosynthetic electron transport 
chain (PETC) by high irradiance may generate excess ROS, 
which can damage various cellular components and result 
in oxidative stress (Apostol et al, 1989; Asada, 1994, 1999; 
Foyer and Mullineaux, 1994). Previously, Han et al. (2012) 
demonstrated that the capability to prevent ROS accumu- 
lation effectively was critical for the survival of H. pluvialis 
under high light stress conditions. 
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Table 3. Isotigs of photosystem- and light-harvesting complex-related genes in the H. pluvialis transcriptome 



Isotig number 



Description 



Accession 



Read ratio 
(red/green) 



isotig04372-04373, 04376-04377 

isotig21 425-21 426 

isotig09349-09350 

isotig26971 

isotig22073-22074 

isotig 16438-1 6439 

isotig 1 2347-1 2350 

isotig2291 5-2291 6 

isotig28347 

isotig09222, 09224 

isotig22287-22288 

isotig39653 

isotig30606 

isotig04580-04581 , 04583-04586, 

04589-04590 

isotig02670 

isotig45545 

isotig35094 

isotig29750 

isotig02655-02657, 02659, 02661 , 
02668 

isotig20097-20098 
isotig31218 
isotig09395-09397 
isotig26863 

isotig03001 -03002, 03006-03007, 
03011-03012, 03016-03017, 
03021-03022, 03026-03027 
isotig041 98-04209 

isotig13851-13852 

isotig 1 9573-1 9574 

isotig05175, 05178, 05180 

isotigOl 744-01 788 

isotig051 76-051 77, 05181-05186 



Chloroplast photosystem I reaction centre subunit III 
(Chlamydomonas sp. ICE-L) 
Photosystem I reaction centre, subunit VIII 
(Chlamydomonas reinhardtii) 
Photosystem I reaction centre, subunit VIII 
(Chlamydomonas reinhardtii) 
Photosystem I reaction centre subunit XI 
(Chlamydomonas reinhardtii) 
Photosystem I subunit O (Chlamydomonas 
reinhardtii) 

Photosystem I reaction centre subunit V, chloroplast 

precursor (Volvox carter! f. nagariensis) 

Photosystem I reaction centre subunit VI, chloroplast 

precursor (Volvox carter! f. nagariensis) 

Photosystem I reaction center subunit psaK, 

chloroplast precursor (Volvox carteri f. nagariensis) 

Photosystem I 8.1 kDa reaction centre subunit IV 

(Chlamydomonas reinhardtii) 

Photosystem I 8.1 kDa reaction centre subunit IV 

(Chlamydomonas reinhardtii) 

Photosystem I 8.1 kDa reaction centre subunit IV 

(Chlamydomonas reinhardtii) 

Photosystem I 8.1 kDa reaction centre subunit IV 

(Chlamydomonas reinhardtii) 

Photosystem I 8.1 kDa reaction centre subunit IV 

(Chlamydomonas reinhardtii) 

10 kDa photosystem II polypeptide (Chlamydomonas 

reinhardtii) 

Photosystem II 47 kDa protein (Dunaliella salina) 
Photosystem II protein I (Dunaliella salina) 
Photosystem II protein V (Dunaliella salina) 
Photosystem II protein Z (Dunaliella salina) 
Photosystem II subunit B (Dunaliella tertiolecta) 

Photosystem II subunit W, chloroplast precursor 

(Volvox carteri f . nagariensis) 

Photosystem II subunit 28 (Chlamydomonas 

reinhardtii) 

Oxygen-evolving enhancer protein 3 
(Chlamydomonas reinhardtii) 
Photosystem II protein D2 (Chlamydomonas 
reinhardtii) 

Major light-harvesting chlorophyll alb protein 2.1 
(Dunaliella salina) 

Major light-harvesting chlorophyll alb protein 2.1 
(Dunaliella salina) 

Light-harvesting chlorophyll a/fo-binding protein 
Lhcb5 (Chlamydomonas incerta) 
Light-harvesting chlorophyll a/b-binding protein 
Lhcb5 (Chlamydomonas incerta) 
Light-harvesting chlorophyll a/b-binding protein 
Lhca6 (Chlamydomonas incerta) 
Light-harvesting chlorophyll a/fo-binding protein of 
photosystem I (Volvox carteri f. nagariensis) 
Light-harvesting protein of photosystem I (Volvox 
carteri f . nagariensis) 



ADQ00182.1 

XP_001 703367.1 

XP_001 703367.1 

XP_001 691 084.1 

XP_001 7001 09.1 

XP_002949489.1 

XP_002956888.1 

XP_002946059.1 

XP_001 70261 1.1 

XP_001 70261 1.1 

XP_001 70261 1.1 

XP_001 70261 1.1 

XP_001 70261 1.1 

XP_001 696588.1 

YP_005089828.1 
YP_005089795.1 
YP_005089824.1 
YP_005089815.1 
AAX76829.1 

XP_002955022.1 

XP_001 690537.1 

XP_001 701 331.1 

1209190A 

ABD91 646.1 

ABD91 646.1 

ABD37916.1 

ABD37916.1 

ABD37908.1 

XP_002957031.1 

XP_002958611.1 



0.59 

0.52 

0.60 

0.63 

0.69 

0.48 

0.52 

0.50 

0.55 

0.49 

0.65 

0.41 

0.69 

0.67 

0.94 
0.63 
0.95 
0.95 
0.96 

0.46 

0.72 

0.77 

0.73 

0.38 

0.31 
0.50 
0.50 
0.44 
0.66 
0.43 
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Table 3. Continued 



Isotig number 



Description 



Accession 



Read ratio 
(red/green) 



isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 
isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 

isoti 



g1 0497-1 0500 

g1 1 1 50-1 1 152 

g23593 

g28923 

g13760 

g13759, 13761 

g20625 

g40410 

g22267-22268 

g29810 

g03544,03547 

g18101-18102 

g03028-03029 

g21 535-21 536 

g06502-06507 

g39761 
g13853 

g13130-13132 

g1 91 49-1 91 50 

g201 07-201 08 

g27041 

g27610 

g27685 

g28381 

g29135 

g29172 

g29506 

g35073 

g0931 3-0931 6 



Light-harvesting protein of photosystem I (Volvox 
carteri f . nagariensis) 

Light-harvesting protein of photosystem I [Volvox 
carteri f . nagariensis) 

Light-harvesting protein of photosystem I (Volvox 
carteri f . nagariensis) 

Light-harvesting protein of photosystem I (Volvox 

carteri f. nagariensis) 

Light-harvesting protein of photosystem 

I (Chlamydomonas reinhardtii) 
Light-harvesting complex A protein (Volvox carteri 
f. nagariensis) 

Light-harvesting complex A (Volvox carteri 
f. nagariensis) 

Chloroplast light-harvesting chlorophyll-a/fo-binding 

protein (Chlamydomonas incerta) 

Early light-inducible protein (Chlamydomonas 

reinhardtii) 

Early light-inducible protein (Chlamydomonas 
reinhardtii) 

Low molecular mass early light-induced protein 

(Volvox carteri f . nagariensis) 

Light-dependent protochlorophyllide reductase 

(Chlamydomonas reinhardtii) 

Chlorophyll a/6-binding protein (Chlamydomonas 

moewusii) 

Chlorophyll a/t)-binding protein (Chlamydomonas 
reinhardtii) 

Chlorophyll a/t)-binding protein of photosystem II 
(Volvox carteri f . nagariensis) 
Major chlorophyll-binding protein (Dunaliella salina) 
Minor chlorophyll a/b-binding protein of photosystem 

II (Chlamydomonas reinhardtii) 
Oxygen-evolving enhancer protein 1 of photosystem 
II (Chlamydomonas reinhardtii) 
Oxygen-evolving enhancer protein 2 (Volvox carteri 

f. nagariensis) 

Oxygen-evolving enhancer protein 2 (Volvox carteri 
f. nagariensis) 

Oxygen-evolving enhancer protein 1 of photosystem 

II (Chlamydomonas reinhardtii) 

Oxygen-evolving enhancer protein 1 of photosystem 

II (Chlamydomonas reinhardtii) 

Chloroplast oxygen-evolving complex/thylakoid lume- 

nal 25.6 kDa protein (Chlamydomonas incerta) 

Oxygen-evolving enhancer protein 2 (Volvox carteri 

f. nagariensis) 

Chloroplast oxygen-evolving protein 3 
(Chlamydomonas incerta) 
Chloroplast oxygen-evolving protein 3 
(Chlamydomonas incerta) 

Oxygen-evolving enhancer protein 2 of photosystem 
II (Chlamydomonas reinhardtii) 
Oxygen-evolving enhancer protein 2 (Volvox carteri 
f. nagariensis) 

Plastocyanin (Scenedesmus obliquus) 



XP_002957416.1 

XP.002948151.1 

XP_002958611.1 

XP_002958611.1 

XP_001 692548.1 

XP_002950368.1 

XP_002947520.1 

ABA01 131.1 

XP_001 695978.1 

XP_001 694751.1 

XP_002958870.1 

XP_001 689464.1 

P22686.1 

XP_001 699932.1 

XP_002951 424.1 

P20865.1 
XP_001 695927.1 

XP_001 694699.1 

XP_002956365.1 

XP_002956365.1 

XP_001 694699.1 

XP_001 694699.1 

ABA01 138.1 

XP_002956365.1 

ABA01 140.1 

ABA01 140.1 

XP_001 6941 26.1 

XP_002956365.1 

P26956.2 



0.51 
0.68 
0.24 
0.29 
0.60 
0.60 
0.03 
0.19 
0.87 
0.71 
0.52 
0.47 
0.39 
0.88 
0.56 

0.10 

0.66 

0.58 
0.67 
0.40 
0.54 
0.83 
1.41 
0.73 
0.76 
0.98 
1.03 
0.00 
0.59 



4332 I Gwakefa/. 



Cells can fight against oxidative stress via several defensive 
mechanisms, including increasing ROS scavenging activity 
or decreasing ROS production. Scavenging of ROS can be 
mediated by a variety of enzymes, such as superoxide dis- 
mutase (SOD), catalase (CAT), ascorbate peroxidase (APX), 
glutathione reductase (GR), monodehydroascorbate reduc- 
tase (MDHAR), dehydroascorbate reductase (DHAR), and 
glutathione peroxidase (GPX). ROS accumulation can be 
also counteracted by non-enzymatic, low molecular weight 
metabolites including ascorbic acid, glutathione (GSH), 
thioredoxin (TXN), and secondary carotenoids (for reviews, 
see Blokhina et al, 2003; Moulin et al, 2010). It is possible 
that the two kinds of antioxidant systems work in concert to 
protect cells from oxidative damage. 

The NGS analysis carried out here identified multiple 
genes coding for GPX (two), GR (three), and APX (two), and 
at least one gene for each enzyme was up-regulated by high 
irradiance, up to ~3-fold (Table 2). As for low molecular anti- 
oxidant proteins, only TXN genes were detected in the NGS. 
A single gene was found for DHAR, the transcripts of which 
were detected only in haematocysts. Two genes identified 
for APX GPX were expressed at a higher level (~ 1.8- to 
3-fold) in haematocysts. All three GR genes were up-regulated 
~1.6- to 3.8-fold by high irradiation. Four genes were identi- 
fied for SOD, among which one transcript displayed a slight 
increase (~ 1.7-fold) in response to high irradiation whereas 
the other three genes were not affected. Likewise, the tran- 
script levels of a single catalase and MDHAR gene were not 
significantly affected by high light conditions. Based on the 
present results, it would seem that an increase in the abun- 
dance of those transcripts related to ROS scavenging activ- 
ity could contribute to the protection of cells from oxidative 
stress induced by high light levels. 

The present NGS analysis revealed that high light strongly 
down-regulated the expression of genes encoding proteins 
involved in light-harvesting activities, such as the light-har- 
vesting chlorophyll a/6-binding proteins of photosystems 
(PS) I and II. High irradiance also decreased the expression 
of genes coding for PSI and PSII reaction centre proteins 
(Table 3). The cytochrome b(f complex functions as a plas- 
toquinine (PQ)-plastocyanin-oxidoreductase to couple PSII 
and PSI, and is also involved in cyclic electron transport. In 
the NGS analysis, seven genes that encode the cytochrome 
b(f complex were identified and four were down-regulated 
under high light conditions. Most notably, expression of the 
most highly active gene in macrozooid cells declined by >30% 
under high irradiance. The transcript level of an identified 
plastocyanin gene also decreased by 55% under high light. 
Among 11 genes encoding oxygen-evolving complexes, the 
three most highly expressed genes were down-regulated up 
to ~2-fold in haematocysts, suggesting that high light condi- 
tions decreased the primary electron-donating activity from 
oxygen-evolving proteins. In addition to its role as a cofactor 
to couple PETC to astaxanthin synthesis, PTOX is known to 
facilitate reoxidation of the PQ pool in the photosynthetic 
electron transport chain independently of cytochrome b(f. 
PTOX can serve as an alternative electron sink, transferring 
electrons from the PQ pool to 0 2 to form H 2 0 (Cournac 



et al, 2002; Joet et al, 2002; Peltier and Cournac, 2002). As 
mentioned earlier, haematocysts displayed higher expression 
of PTOX genes. All of these transcriptomic differentiations 
that serve to reduce the light-harvesting, reaction centre, 
and linear electron transport activities may help to relax 
over-reduced PETC and lead to a decrease in ROS forma- 
tion. The amounts of ROS produced from late exponential 
Haematococcus macrozooid cells were monitored following 
exposure to high light conditions for 4 d using H 2 DCFDA. 
Consistent with the transcriptomic adjustment which was 
observed, the production of ROS decreased dramatically in 
cells after 3 d under high irradiance and remained at lower 
levels thereafter (Supplementary Fig. S8 available at JXB 
online). 

In summary, the transcriptomic and lipidomic analyses sug- 
gested that under high irradiance stress, H. pluvialis cells were 
protected from photooxidative stress by multiple defensive 
systems. High irradiance induced formation of astaxanthin 
and TAG, while reducing membrane glycerolipids, espe- 
cially those glycolipids making up the photosynthetic com- 
plexes and chloroplast membrane matrix. The transcripts of 
genes involved in astaxanthin biosynthesis — including those 
responsible for supplying the isoprene precursor via the MEP 
pathway or converting IPP to astaxanthin — were globally 
up-regulated to a comparable extent. Accumulation of TAG 
under high irradiance was attributable to the remarkable 
up-regulation of de novo fatty acid biosynthesis at the gene 
expression level as well as to the moderate elevation of the 
TAG assembly pathways. Consequently, enhanced astaxan- 
thin and fatty acid biosynthesis consumed excessive electrons 
from over-reduced PETC under stress. Increasing enzyme- 
based ROS scavenging activity and reducing light-harvesting 
capacity, primary photochemical energy production, and lin- 
ear electron transport (as inferred from the transcriptomic 
analyses) also augmented the ability of H. pluvialis cells to 
survive under photooxidative stress. 

Supplementary data 

Supplementary data are available at JXB online. 

Figure SI. Minor and trace triacylglycerol (TAG) species 
in H. pluvialis. 

Figure S2. Profile and quantification of PG (A), PC (B), 
DGTS (C), PI (D), and PE (E) in H. pluvialis. 

Figure S3. The length and coverage of contigs, isotigs, and 
reads from the H. pluvialis transcriptome database. 

Figure S4. A detailed analysis of GO terms of the three 
functional classes of H. pluvialis. 

Figure S5. Correlation plot of housekeeping genes in 
H. pluvialis transcriptomes. 

Figure S6. KEGG pathway analysis of H. pluvialis show- 
ing differential gene expression (>2-fold) between macrozoo- 
ids (A) and cysts (B). 

Figure S7. Mevalonate (MVA) and 2-C-methyl-D-erythritol 
4-phosphate (MEP) pathways in H. pluvialis. 

Figure S8. ROS production of H. pluvialis. 
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Table SI. Abbreviations used in this paper together with 
their full names. 

Table S2. Total carotenoid and chlorophyll content under 
the different culture conditions. 

Table S3. Housekeeping genes in the H. pluvialis 
transcriptome. 

Table S4. Membrane glycerolipid biosynthesis-related 
genes in the H. pluvialis transcriptome. 
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